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Abstract 

We scrutinized hadronic Regge trajectories in a framework of two different models 
- string and potential. Our results are compared with broad spectrum of existing 
theoretical quark models and all experimental data from PDG98. It was recognized 
that Regge trajectories for mesons and baryons are not straight and parallel lines in 
general in the current resonance region both experimentally and theoretically, but very 
often have appreciable curvature, which is flavor-dependent. For a set of baryon Regge 
trajectories this fact is well described in the considered potential model. The standard 
string models predict linear trajectories at high angular momenta J with some form of 
nonlinearity at low J. 



Introduction 

In the past few years Regge trajectories (RT) have attracted the interest of many authors, 
who build quark models of baryons, mesons, glueballs and hybrids [1-19]. Quite a number 
of approaches have been used to attack this problem: the WKB-approach |l|, |2|] , the Wilson 
loop model g |, the ^-expansion technique [f|, q-deformed algebra approach ||, 1/N C - 
approach [7|, spectrum generating algebra (SGA) model ||, Filipponi model ||, NRQM 
jTOfl , the string models fllf - Hl5[] , extended covariant oscillator quark model for glueballs 



|T§, iV/D-method [H]] and others. 

Probably, the most interesting questions, which are under investigation, are the follow- 
ing: (a) Are the RT really straight lines in the entire energy interval, or this is only true 
asymptotically? (b) Do the trajectories for mesons, baryons, glueballs and hybrids have the 
same slope? (c) What is the flavor-dependence of RT? (d) What is the intrinsic connection 
between kinematics, the type of the potential and the straight RT? (e) Is it possible to 
connect the linear RT with parity-doubling phenomena (tower structure)? (f) What is the 
effect of color screening on heavy and light quarkonia and baryon's RT? (g) When does the 
asymptotic regime (J=?) really start for baryons, mesons, glueballs, hybrids for both parent 
and daughter (ancestor) trajectories? (h) What is the dependence of the character of RT on 
the scalar/vector structure of the confinement potential? 

Different groups pursued variety of models and different aspects of this problem. It has 
been established long time ago, that the experimental RT for N, A baryons are not strictly 
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straight lines fI8 |. Authors of this seminal review and Hendry in review [19j considered the 
facts of nonlinear behavior of the RT in mass- squared. Rather, as Hendry concludes, baryon 
resonances seem more linear as a function of cm momentum. 

In this paper we will show that both experimentally and theoretically RT may deviate 
from straight and parallel lines, depending on particular family of baryons, mesons, glueballs 
or hybrids and energy region. We will scrutinize this problem in a framework of two different 
models — string |0|, and potential [JTD | . Then we will contrast our results with broad 
spectrum of existing theoretical models and experimental data for different flavors and draw 
the conclusions. 

1. Regge Trajectories in the Potential Quark Model 

In our recent series of papers |10[, based on hamiltonian (1) and the method of hyper- 
spherical functions (HF) 
was given. 



(J, a description of N, A, fl resonance spectra and partial widths 
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we introduce the constants a s , C a and C t , which determine the 
strength of the Coulomb, contact and tensor potentials respectively. The use of the hamilto- 
nian (1) allows us to obtain better agreement with experiment for resonances of positive and 
negative parity, and also to describe resonances with both large J and M, and with small J 
and M. 



We showed |10| that it was appropriate for such a description to take advantage of the 
concept of yrast states and yrast lines from the theory of atomic nuclei rotational spectra, as 
well as to make use of the concept of RT. One of the main results was that both theoretical 
and experimental spectra are nonlinear trajectories in Chew-Frautschi plots throughout the 
experimental region. Our model predicts a whole series of high-lying N, A resonances, 
represented in Table 1. 

Table 1. Masses of the predicted N, A resonances in the model [M 



N 11/2+ (2.39) 


N 13/2- (2.90) 


N 15/2+ (3.00) 


N 17/2- (3.45) 


N 19/2+ (3.45) 


N 21/2- (3.95) 


A 15/2" (3.48) 


A 17/2+ (3.49) 


A 19/2- (4.00) 



The present paper generalize the model [jnj to u, d, s, c, b flavors (it is impossible to 
create top baryons and mesons) and a wider range of angular momenta L = -r- 20. We 
obtain the spectra of BR in the same way as in [10||, by solving Schrodinger equation (SE) 
with the HF method. 

When the hadron wave function (WF) is expanded in the HF basis and substituted in 
the SE, one generally finds an infinite system of differential equations for the radial WF 



(RWF). However, as was shown in [20] for a system of identical u,d-quarks, the coupling of 
channels is weak, therefore it is sufficient to include only several terms in K, grand orbital 
momentum, in the WF expansion. With increasing excitation energy, the contribution of 
Hhyp vanishes, hence the coupling of channels in general can be neglected. On the other 
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u, s, c, 



b), 



hand, as we will consider only symmetric baryons with increasing quark mass 
the above mentioned arguments will apply even more strongly. 

Now, let us introduce Jacobi and hyperspherical coordinates for a three-body problem 
The Jacobi coordinates are defined as usual (please see [|Hj] for all the details) 



R=(r 1+ r 2 + r 3 )/3, r) = (n - r 2 )/ v / 2, £ = (^/3)[(r! + r 2 )/2 - r 3 ]. (2) 

The hyperspherical angle 9 and radius p are defined as follows 

n = pcos6, £ = psinfl, < 6 < vr/2, p = (p 2 + £ 2 ) 1/2 , < p < oo. (3) 

So, our WF and q-q potentials will depend on new variables, r] and £. 

We will work in the hypercentral approximation (HCA), where only that part of the inter- 
action which is invariant under rotation in six-dimensional space (p, f2 5 ) is taken into account. 
Because in HCA baryons are pure rotational states in six- dimensional space, and because a 
very soft potential is used in our model, no strong correlations generating a so-called quark- 
diquark state can occur. So, we will neglect the 6*-dependence of our q-q potential. The 
introduction of ^-dependence in the potential mixes the states belonging to one value of L, 
and slightly shifts their positions in the baryonic mass spectrum. As was proved by Richard 
the most commonly used potentials in hadron spectroscopy (including ours) are very 



20 



close to hypercentral (please see |]TU[ for all the details). 

2. Mass Spectra and RT in the Potential Quark Model 

The method of 
but we will briefly recapitulate 



Now let us discuss our results for mass spectra and RT of baryons. 
numerical solution of the SE was described in detail in 
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it here. The system of differential equations (SDE) was reduced to a system of first -order 
differential equations, which we integrated by determining a full set of Cauchy solutions. We 
then constructed the required solution by imposing the following boundary conditions on 
the RWF: Fj(p = 0) = and Fj(p = R^) = 0. The eigenvalues (EV) were determined by 
zeroing the determinant constructed from the fundamental system of solutions for SDE. For 
the parameter we always choose a much larger value than the radius of the corresponding 
excited state (p)j, which grows with L and N r (the radial quantum number). We always 
choose Roo such that if we take R^ = 2/?^, then the corresponding EV of any excited state 
will not change more than 1%. We use in our computation the following set of parameters 
(Table 2). With this input we have calculated the mass spectra, RT, mean radii (p)j and 
slopes for u, d, s, c, b flavors, N r = 0, 1, 2 and momenta range L = 4- 20. 

Table 2. Input parameters in the model pi 



a s = C a = 1 


C t = 0.2 


b = 350 MeV/fm 


V = -513 MeV 


m u = 330 MeV 


m d = 330 MeV 


m s = 607 MeV 


m c = 1500 MeV 


m b = 5170 MeV 



We note that the solution of the SE in the single-channel approximation with the cen- 
trifugal potential V c = (K + 3/2) (K + 5/2) /(2mp 2 ) is equivalent from the mathematical 
point of view to the Davydov-Chaban model p2"] , in which the rotational spectrum of a 
non-spherical nucleus with variable moment of inertia is calculated. Because the moment of 
inertia is not constant but varies linearly with the angular momentum J of the nucleus, the 
rotational spectrum, instead of the quadratic behavior J(J+ 1), it would have for a constant 
moment of inertia, approaches one linear in J. 
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Approximately the same thing happens with the hadron: for large excitation energies the 
spectrum becomes linear in J, and the hadron becomes a strongly extended system along 
its symmetry axis. This should lead to a linear nature of the baryonic mass spectra along 
yrast lines. 



Our predictions are in accord with findings by Hey jl8| and Hendry 
Let us start a detailed, sector by sector comparison of our results for different flavors. 
Note that RT (M 2 = M 2 ( J)) for baryons have different types of curvatures: some are convex; 
others are concave functions of J (see Figs. 1,2). The u-d (N, A) trajectories are concave for 
N r = 0, 1, 2 (Fig. 2), whereas strange (sss) RT for N r = 0, 1, 2 are "stereotypical" — they 
are nearly straight lines with slowly varying slopes a' (Fig. 3). Charmed (ccc) and bottom 
(bbb) trajectories are convex for N r = 0, 1, 2. 

When we analyze the mass spectra M = M(J) for u-d, s, c, b-baryons, they are all 
convex functions of J, with a rather complex dependence on m q (flavor) and N r . 

If we fix m q and look at the Chew-Frautschi plot for N r = 0, 1, 2, we can see that these 
three curves generally are non-parallel (Fig. 2), that strongly contradicts the conventional 
picture of the parallel RT. It is noteworthy that Olsson et. al. [[HJ recently examined the 
meson sector, using a relativistic flux-tube model, and noticed nonlinearity of RT at low 
angular momenta J. But it remains unclear whether their model accounts adequately for 
physical observables, because it describes a very limited number of states in the mesonic 
sector. 

Slopes for the u-d family start above 1 GeV™ 2 , with the largest value for N r = 2; the 
curves decrease almost monotonically with changing the rate of decreasing near L — 4, for 
N r = 0,1,2 (see Fig. 3). 

Slopes for the s-family differ from all the other cases, because they are rather weak 
functions of L. Slopes for daughters N r = 1, 2 start just above unity, then slowly decrease to 
0.87 GeV -2 and 0.78 GeV -2 , respectively, whereas the parent slope fluctuates slightly over 
the interval, spanning the range 0.97-0.91 GeV -2 (see Fig. 3). 

Slopes for the charm family are highly nonlinear functions of L. The N r = 0, 1, slopes 
increase monotonically, starting from 0.39 GeV -2 and 0.55 GeV -2 and approaching 0.889 
GeV -2 and 0.893 GeV -2 , respectively. The N r = 2 slope has a dip at L = 4, and grows 
continuously, reaching 0.897 GeV -2 . These values are almost identical to Olsson's results for 



mesons for asymptotic L |13| . 

Bottom baryon slopes differ sharply from the u-d-s-c sector, first, by their small magni- 
tude, and second, by the significant increase of the slopes along the trajectory (the B-slopes 
increase by about an order of magnitude, from 0.055 GeV -2 to 0.504 GeV -2 ). Daughter 
slopes for N r = 2 fluctuate, but still increase with L. 

For convenience, we present the set of median values (ct')i for the whole flavor multiplet 
(see, Table 3). It is interesting to note that the median values of a[ are almost independent 
on N r for the u-d and s-families. 

Table 3. Median values (a')i for trajectories with N r = 0, 1, 2 



N r 


Up-Down 


Strange 


Charm 


Bottom 





0.90 


0.95 


0.73 


0.31 


1 


0.87 


0.95 


0.77 


0.34 


2 


0.84 


0.94 


0.81 


0.37 



The expectation values of the hyperradius (p)i are basically smooth increasing functions 
of L and N r , and decreasing functions of m q . 

We proved that the slope of the trajectories decreases with increasing quark mass in the 
mass region of the lowest excitations. This is due to the contribution of the color Coulomb 
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interaction, which increases with mass, leading to a curvature of the trajectory near the 
ground state. In the asymptotic regime, the trajectories for all flavors are linear and have 
the same slope a' ~ 0.9 GeV~ 2 . 

After careful numerical evaluation of baryonic and mesonic spectra for all flavors we have 
shown that the point of establishment of linear RT depends on the exponent v in the power 
law potential in hamiltonian (1), and occurs at larger L with larger v. For an oscillator 
confinement [y = 2) linear regime started only from L > 20, for linear confinement {y = 1) 
linear regime started from L > 18. 

We proved that for mesons linear RT is established earlier than for baryons as a function 
of L. The reason lies in centrifugal energy term, which has LiL + 1) dependence for mesons 
and (L + 3/2) (L + 5/2) dependence for baryons. As we see, three-dimensional corrections to 
L(L + l)-law are important for baryons. 

3. String models 

The string models are widely used for describing orbitally excited hadron states by virtue 
of some remarkable features: (a) the direct analogy between the string with linearly growing 
energy and the QCD confinement mechanism of connecting quarks (antiquarks) by the gluon 
field tube [p^j ; (b) the strings are relativistic by definition^ (c) the energy E = M and the 



angular momentum J of a rotating open (massless) string are connected by the Nambu 



relation P4" 



J = a'M 2 , a' = (27T7)- 1 , (4) 

where 7 is the string tension. This fact allows us to apply the string models to describing 
the RT. 

The massless string generates the strictly linear RT (|4]). But the meson model of rel- 
ativistic string with massive ends |TI[] (more realistic one in comparison with the massless 



case) results in the more complicated behavior of RT. For this model the relation (|4]) takes 



place only in the high energy limit ITT], |12jl . 

Shortly after this meson model string models of baryon were suggested in four variants 
differing from each other in the topology of spatial junction of three massive points (quarks) 
by relativistic strings: (a) the quark-diquark model q-qq [12] (from the point of view of 
classical dynamics it coincides with the mentioned meson model of relativistic string with 
massive ends ]TT], |2B|); (b) the linear configuration q-q-q with quarks connected in series 



15 , 25] ; (c) the "three-string" model or Y-configuration with three strings joined in the 



fourth massless point (junction) [ER E7| and (d) the "triangle" model or A-configuration 



that could be regarded as a closed string carrying three pointlike masses |^8|, p9| . 

All mentioned string hadron models may be constructed on the base of the action JTT 

H> El 

T 2 <Tjv(T) N 

s = ~Jv I V^O'-^^^ + E^V^V)}^- ( 5 ) 

n <7*(t) i=l 

Here rrii are masses and N is the number of the material points (quarks, antiquarks or 
diquarks)^, N = 2 for the meson and quark-diquark models and N = 3 for other models 
of baryon, x M = -X" M (r, a) is the string world surface in rf-dimensional Minkowski space with 

lr The relativistic string dynamics results from the extremization of a world surface area swept by the 
string in Minkowski space. 

2 For brevity we will use the term "quark" instead of "material point" , keeping in mind that the action 
(ra) on the classic level doesn't describe the spin and other quantum numbers of quarks. 
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signature +, 



X" = d T X", X'k = d a X", c = 1, h = 1; V? = j-X*{T,<Ti(T)) is the 



for all 



tangent vector to the i-th quark trajectory a = 0j(r), <J%{t) < a i+ i(r), 
models, except for the A configuration. 

In the "triangle" model the equations a = a* = <j (t) and a = o" 3 (r) define the tra- 
jectory of the same (the 3-rd) quark. This fact may be written as the closure condition 
X^fr, cTo(r)) = X M (r* , cr 3 (r*)) f28fl . The parameters r and r* aren't equal in general. 

The equations of motion and the boundary conditions on the quark trajectories in these 
models are deduced by variation and minimization of action (|5]). In particular, under the 
conditions of orthonormality X 2 + X' 2 = 0, {XX 1 ) = (they may be obtained for all 
configurations [^3|, |28|) the equations of motion become linear 

(6) 



X» - X"^ = 0, 
and the boundary conditions take the simplest form 

-Ft. 



d Vt 



in. 



dr 



(7) 



o'Ar 



for a quark at an end of the string, this expression 



X»\ 



=o-i+0 



7 [X'" + <(t)X"] 



for a quark 



On this 



Here if = ± 1 [X'» + g'^t) X"] 
takes another form Ff 1 = ^[X'^ 

in the "triangle" model and for the middle quark in the q-q-q configuration |f28 
trajectories the derivatives of X M (r, a) are not continuous in general. 

The string models are successfully applied to describing the main or parent RT. For this 
RT, that is for the orbitally excited hadrons the rotational motions of all string configurations 
(flat uniform rotations of the system) are used Ul2 , 15, 







The solution of this type satisfying Eq. (|6]), conditions (0) and describing the uniform 
rotation of the rectilinear string is well known for the meson string model |TT1 , f23f (or its 
equivalent q-qq model) and for the q-q-q baryon configuration. It may be represented as 
(N — 3 for the q-q-q) 



X° = t = r, 



X 1 + iX< 



sin ua ■ e 



(8) 



Here u is the angular velocity, X 1 = x, X 2 = y, a G [a±, <tn], <?i = const, o\ < 0, ajv > 0. 

For the linear q-q-q configuration the middle quark is at rest at a = oi = 0. But this 
motion is unstable with respect to centrifugal moving away of the middle quark that results 
in the complicated quasi-periodical motion with varying distance between the nearest two 
quarks ||25|| . However the minimal value of the mentioned distance does not equal zero, in 
other words, the system q-q-q is not transformed in the quark-diquark one as was supposed 
in Refs. 
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So the q-q-q system is probably applicable not to pure orbital excitations but 
to radial ones. 

The rotational motion of the "three-string" model with the junction at rest and with 



rectilinear string segments joined in a plane of rotation at the angles 120° |26|, |27j is described 
by the expression similar to Eq. (|8|). 

For the meson string model with massive ends and for the baryonic configurations q- 
qq, q-q-q and Y the energy E = M and the angular momentum J (z projection) of the 
considered rotational motions are 



11, 12, 15 



M 



N r 

E 

i=l 



7 



arcsin t>j + 



+ AM, 



(9) 



J 



N 

E 

t=i 



2uj\ 



arcsm Vi 



1 



+ Si 



(10) 



where the velocities of moving quarks = sin \ua, 



+ 1 - 



The presence of 



III , LP 

S) is taken into account as the correction AM 



quark spins with projections s, (X i=1 
to the energy of the classic motion. In Refs. []12], |TjJ] this correction is due to the spin-orbit 
interaction (the spin-spin correction is assumed to be small in comparison with the spin-orbit 
one at high J). 



L. D. Soloviev [14| in the frameworks of another approach constrains the string dynamics 
only rectilinear string motions (H) (with motions of the cm.) and introduces into action (|5|) 
the quark spin terms with some additional terms. The latter ones compensate the influence 
of the quark spins so that the resulting motion is the same as in the spinless case. Such a 
form of compensation results in Ref. |TJj in the same expressions (^|) and (|TU[) with AM = 
and two corrections in Eq. fllPl) . The 1. h. s. of Eq. fllPl) in Ref. |14| is substituted by 
J( J + 1) due to some form of quantization and the spin term S = Yli=x s « ^ s substituted by 



the phenomenological parameter a n resulting from the quark interaction at short distances. 
The values a n (different for various RT and depending on J for heavy quarks) are calculated 
by fitting. 

Using this approach L. D. Soloviev |14j] describes a lot of meson states, in particular, the 
states with zero orbital momentum (tt, p), where the string models are unlikely adequate. 

The rotational motions (uniform rotations about the system center of mass) for the 
"triangle" baryon model jLR [28|, |29| may be presented in the form 



T 



X 1 + iX 2 = u(a) 



Here the values <jj, D = cr 3 — cr , X 2 (r, dj) = V 2 , t* — r = T are the constants, the complex 
function u(a) = cosua + Bi sinua, a G [cj, CTj+i] is continuous in [ao,^]. Owing to 
Eqs. (0) six complex constants Ai, Bi are proportional to A\ (here we put o\ = 0) 

A = A u B = K/i! - iK) A 1} B l = -\{h x + %K) A u 
A 2 = [l + i/i 2 rfi(2ci - Mi - iKdx)] A x , B 2 = -\[2h 2 c\ + (1 - h^cxd^hx - iK)] Ax, 

where |^4i| 2 = (1 — V 2 )/u 2 = 2T/{KDuj 2 ) and we use the following notations: 

K = 2d 2 dQ 1 d{ 1 sin uT/(G 2 G3 — 1), di = sincj(cr i+ i — <7j), q = cosu(ai + i — <Ji), 
Gi = (hidi-idi - di_ x Ci - c^xdi) / d i+l , h = um^-^Vi^ 1 . 

The notations q, d i: hi, Gi are cyclical: c i+3 = q, d i+3 = di . . . 

Expression ( pTT]) is the solution of Eq. (||) and satisfies the orthonormality, closure and 
boundary (|7|) conditions if the parameters are connected by the relations |29|, |30| 



2cos^T = Gx + G 2 + G 3 -GxG 2 G 3 , D/T + T/D = (K 2 + 4 + h\)/{2K), 
(Cj+i — Gi) di = (GiGi + i — l)(dj_ 1 Cj + i — c.;_i<ij +1 ). 

For given m^, 7 and a parameter measuring the rotational speed (for example u) one can 
find all mentioned values by solving these equations. 

A string position for this state (a section t = const of the surface (|Tl~D) is the curve 
composed of three segments of a hypocycloid joined in three points (the quark positions). 
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Hypocycloid is the curve drawing by a point of a circle that is rolling inside another fixed 
circle with larger radius. 



The energy and the angular momentum of this state are p8| , p9 

3 



M 



70(1 



mi 



D 2 



i=l 



1 



V 2 



+ AM, J=±- 
2uj 



rriivf 



D 2 



v? 



+ S. (12) 



29 



where the quark velocities Vi = ^Tdi-idi(Gi-\Gi+i — l)/(Dd i+ ismujT). 

A set of topologically different configurations of the system is classified in Refs. 
It is described, in particular, by integer parameters n and k which (if there are no singularities 
X 2 = in the segment a G [ctq, a\}) are determined by the following way: 



D 



n 



lim 



k = lim 



T 



"li^O O"! — (To 



Here k — n — 2, n — 4, . . . , — n + 2. In the case n = 2, k = ( "quark-diquark" state) two 
quarks coincide, one of di equals 0, T = 0. For this state the world surface (|TTD reduces to 
(|^) and the model "triangle" — to the q-qq one with the tension 27 (rectilinear segment is 
the particular case of the hypocycloid). In the state n = 3, k — 1 (simple state) the string 
has a (curvilinear) triangle form, for other n and k (exotic states) it has more complicated 
form with the cusps moving at the speed of light ||29|, [3(| . 

The simple states of this model are used below for describing the RT. It was shown in 
Ref. [|3T] that they are stable (small disturbances don't grow unlike in the q-q-q case |f2~5|), 
in particular, with respect to transformation into the "q-qq" state with n = 2, k = 0. 

In the two variants of describing RT with the help of string models in Refs. |T^] and JE) 
the spin-orbit correction to the energy 



AM SI = £/3(«i)(w«i 



(13) 



were used in the two different forms. In Ref. (fTS] the expression (3{vi) = —[(1 — v 2 ) x l 2 — 1] is 
due to the Thomas precession of the quark spins. It is obtained under the assumption that 



in the quark rest frame the field is pure chromoelectric [0, |32fl. The alternative assumption 
about pure chromoelectric field in the rotational center rest frame results in (3{vi) = 1 — (1 — 

v 2 y/ 2 MM- 



The ultrarelativistic asymptotic behavior of the dependence J(M) for the q-qq or mesonic 
(m), Y and A configurations has the form 



J ~ a'M 2 



N 

uM 1/2 rri 
i=i 



3/2 



AJ, 



(14) 



where a' = 

2-7T7 

equals A J = 




m, 
Y, 
A, 



-3/4 



m, 
Y. 
A. 



The term A J 



1/2 for 



2/3, 

(2/3) 3 / 2 , 

k 2 ), A, V^7 [ Vrn {n 2_ k 2 ) 

=1 Sj[l — for the spin-orbit correction (|T3|) [ |i~5|l and AJ 

the string model fl4jl . 

The Regge slope a' ~ 0.9 GeV -2 is close for mesons and baryons. So the effective value 
of string tension 7 is to be different for the baryon models Y, "triangle" (simple states) and 
q-qq. The following values of the parameters are used here and in Ref. |15[ (for this values 
of 7 the string baryon configurations q-qq, Y and A result in very close RT): 



s 



Table 4. Input parameters in the string model IT 



7 = 7 9 _ gg = 0.175 GeV 2 


1y = §7 


is = §7 


m u = md = 130 MeV 


m s = 300 MeV 


m c = 1600 MeV 



In Refs 
model the values m U( i 
term in the correction 



TJ| for describing baryon (and meson) RT with the help of the quark-diquark 
= 340 MeV, m s = 440 MeV are used. But the Thomas precession 
OP results in the diquark mass strongly depending on the diquark 



spin state (this mass varies from 220 to 550 MeV). On the other hand in Ref. [TJ| the current 
quark masses are used for describing the mesonic RT. 

In this paper we use the string models (q-qq, Y and A for baryons) with the parameters 
in Table 4 and the spin-orbit correction ([Hj) with (3{vi) = 1 — (1 — vf) 1 ^ 2 for describing parent 
RT. 

The results for light baryons (N and A RT) are shown in Figs. 4 and 5. The particle 

In Fig. 4a the N-baryon states with J p 



data are taken from the PDG98 issue |35 
3/2 



1/2 H 



5/2 + , 7/2~ . . . are placed and described by our two models as two different RT. The 
red dotted line shows the potential model (TTJ predictions for the RT with positive parity 
p l/2 + , iV(1680) 5/2 + . . .; blue dash-dotted line corresponds to the blue marked states with 
negative parity. The same notations are used in Figs. 5-8. 

The dependence J = J(M 2 ) for the string model "triangle" and "three-string" is shown 
in Figs. 4-7 as black solid lines (under conditions in Table 4 these curves for the both 
configurations practically coincide). Dots correspond to the quark-diquark configuration. 
The results of the latter model are very close to that of two previous ones. 

String hadron models [fL2]| — |15| at the modern stage are (semi) classical and do not de- 
scribe the isospin and parity. So in this paper we use only the spin projection 5 = s i (here 
J = L + 5) for modeling various RT. For example, in Fig. 4a the N states with positive and 
negative parity may be described as two different RTQ in all (potential and string) models. 
In the string models the spin states 5 = 1/2 and 5 = —1/2 are used for these two RT - 
the curves J = J(M 2 ) are close to rectilinear and fit the majority of states, except for the 
nucleon (L = 0) and the states with M > 3000 MeV which need confirmation (they are 



omitted from the summary table p5|). But the latter baryons are described in the potential 
model |nj rather well. 

For the A-resonances in Fig. 5 the similar picture take place. The states on the A(1232) 
trajectory with positive parity are supposed to have 5 = 3/2. 

In Fig. 4b the N-baryon states with J p = l/2~, 3/2 + , 5/2~ ... are described by two 
different RT in the potential model |10| or by one RT with 5 = 1/2 in the string models. 
Remind that the string models are not applicable to the states with small J or L. Otherwise, 
the potential model (TJJ is the most adequate for small L. It describes the strongly nonlinear 
RT generated by N(1535) l/2~ (Fig. 4b) and A(1910) 1/2+ (Fig. 5b). The string models 
don't predict such nonlinear behavior (convexity) at small J. 

In Fig. 6 the results of the string models A and Y (q-qq is close to them) for the strange 
baryons A, E and H are represented. The potential model |TJJ predictions are calculated 
for Q baryons. The description of the string models is the most adequate for the particles, 
whose spin state may be interpreted as 5 = ±1/2. For the RT with 5 = 3/2, for example, 
E(1385), Q(1670), A(1232) in Fig. 5 the curves lie higher the lowest states with L = and 

3 One can describe all states in Fig. 4a as one RT jl2j with the help of various string models. But such 
an approach requires some (depending on parity and isospin) corrections to energy and is not applicable to 
other trajectories in Figs. 4b, 5. 
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L = 2. For larger J = L + 3/2 the correspondence is better (Fig. 5) but we have no reliable 
strange E, H and f2 mesons with J > 11/2. There are some reasons of such results in the 
string models ||15|| : (a) they are not adequate for low L, (b) the spin-spin interaction must 
be taken into account at low L, (c) perhaps, the form of the spin-orbit correction (|i3|) needs 



some improvement. 

In the charmed sector (Fig. 7) we know only a few states with J = 3/2 and only one 
of them — A^(2625) - - may be interpreted as the orbital excitation of A+. In Fig. 7 
these particles are described by the q-qq, Y and A string configurations with the spin state 
S = 1/2 and the summary spin of the light quarks s u + Sd = (the upper curves in Fig. 7) 
similar to the strange A RT in Fig. 6. For the state E c the trajectories with S = —1/2 and 
s u + Sd = —1 is suggested. For other charmed baryons with L = the considered string 
models are unlikely adequate. 

For any string baryon model the corresponding meson-like model with the same (or 
corresponding) string tension, effective quark masses, types of spin or other corrections must 
describe the RT in the meson sector. The results of this "meson test" for the string model 
are represented in Fig. 8 for light and strange mesons in comparison with the results 



of the string meson model [|14j]. The latter model contains more fitting parameters (a set of 
a n in addition to the common ones 7 and rrii) so it describes some mesons better. But the 
results of the model 13] are satisfactory for the majority of light and strange meson states. 



Conclusion 

During the last two decades, a number of authors have shed some light on the problem of 
RT in the meson, baryon and glueball sectors. Our approach in this paper is to analyze both 
pure experimental RT and predict properties of the RT in two different models: string and 
potential. This way we hope to understand better the fundamental property of RT — its 
linearity and when it is broken. In particular, our analysis of all available data from PDG98 
55| reveals that the following RT are essentially nonlinear (see Table 5). 

Table 5. Slopes for nonlinear baryon, meson RT 
(a', average (a'), mean square deviation o in GeV~2) 



RT for baryons 


slopes a' for neighbor pairs 


(a') 


a = ^D 


N 3/2" parent 


0.80 1.02 0.702 0.45 


0.74 


0.20 


A 1/2+ parent 


3.53 1.20 0.40 


1.71 


1.33 


N 1 /2~ parent 


4.44 0.89 


2.67 


1.78 


N 1/2- radial 


2.73 0.61 


1.67 


1.06 


N 3/2" radial 


1.72 0.70 


1.21 


0.51 


A 1/2" parent 


1.82 0.98 1.11 0.60 0.34 


0.97 


0.51 


RT for mesons 


slopes a' 


(a') 


a 


f(0 ++ ) parent 


3.00 0.78 0.94 


1.58 


1.01 


K(CT) parent 


0.73 1.14 0.34 1.25 


0.87 


0.36 


p — a 2 parent 


0.87 0.90 0.82 0.69 2.08 


1.07 


0.51 


T radial 


0.09 0.15 0.21 0.16 0.30 


0.18 


0.070 


Xb(lP) parent 


1.58 2.37 


1.97 


0.40 


X&(2P) parent 


2.02 3.66 


2.84 


0.82 


J/\I> radial 


0.25 1.60 2.11 1.03 0.46 


1.09 


0.69 


fo radial 


1.09 2.68 0.55 6.1 1.68 


2.42 


1.97 


f2 radial 


2.39 2.48 4.27 1.64 1.82 4.95 1.9 1.42 5.13 


2.88 


1.39 


X C (1P) parent 


1.54 3.15 


2.34 


0.80 



in 



So, in total we have 6 baryon and 10 meson experimental RT, which we consider as 
essentially nonlinear. We definitely witness the fundamental fact, that RT are not linear as 
a rule — it depends on intrinsic quark-gluon dynamics. One can observe nonlinearity for 
various RT in all current resonance region — for low or large momentaP]. Slopes for many 
experimental RT strongly differ from the standard a' ~ 0.9 GeV~ 2 and vary significantly 
along the given RT (a ~ 1 GeV~ 2 ). 

Now we return to theoretical predictions and interpretations of RT for baryons, mesons 
and glueballs. In the series of papers ||, the authors introduced a new procedure for the 
solution of the SE for mesons, which is based on the expansion in the Planck constant h. 
In contrast to a number of papers, which employed a WKB approximation, the authors of 
|| assumed that their procedure will work for large as well as for small values of the radial 
quantum number N r , and that this method is complementary to the WKB method. The 
authors investigated p — a 2 , J/\l/ and T parent and daughter RT and show that all usually 
employed potentials lead to nonlinear RT in the resonance region. However, authors [0] did 
not consider baryonic systems and they did not show that the model adequately describes 
the experimental meson spectrum, which diminished the value of their papers |>J. 

The authors of investigated the baryon's RT in a relativistic approach, based on the 
method of vacuum correlators (MVC). While deriving the solutions of the dynamic equations, 
the authors neglected the spin, isospin, Pauli principle and one-gluon exchange potential, 
and thereby deduced that the RT were linear and found the N, A mass spectrum up to 
and including K — 6. They found the slope of baryonic RT to be equal to the slope of 
mesonic RT (a' = 0.75 GeV~ 2 ). It is interesting that authors compare the results of this 
approach with the previous paper of one of the authors (F.R.) 0, which was based on the 
nonrelativistic SE with the power law confinement potential r 2//3 . The authors of || came 
to the conclusion that both methods gave very close spectra and linear RT. However, the 
above-cited authors |J sharply criticized the paper for misrepresenting the real picture 
in the resonance region, where the trajectories have a nonlinear character. Simonov, the 
leader of method 0, lends support to our thesis in his collaboration with new authors; 
they very recently have noticed a weak nonlinearity of the mesonic RT |4[] using the same 
approach (MVC). Their observation of the nonlinearity intersects that of Olsson et al [ |T3|j , 
who analyzed the mesonic RT in the relativistic flux tube model. 

J.Dey et al |J described RT for mesons and baryons in g-deformed model and found 
a strong effect of nonlinearity for these RT, for rotational and radial excitations, including 
daughters. They study tt — b, p — a, K, K*, (p, u — f, T, \& meson families and N, A, A, £ 
baryon families, computing spectra up to J = 17/2. 

In his seminal paper fTj 't Hooft developed an 1/N C relativistic, toy model for mesons, 
and he got a RT which possess some nonlinearity: "... deviations from the straight line are 
expected near the origin as a consequence of the finiteness of the region of integration and 
the contribution of the mass terms..." Mostly, the results of numerical analysis are described 
by the following WKB mass formula with nonlinear logarithmic term 

M 2 (n) vr 2 n + («i + a 2 ) logn + C st (a 1 ,a 2 ). (15) 

Iachello et al || have noticed nonlinearity in the experimental masses for T families 
"...We next plot in Fig. 13 and 14, the "vibrational" RT for the \l/ and T families. We note 
that experimental masses no longer fall exactly linear on trajectories, but their trajectories 
are slightly bent, a feature also anticipated in the 't Hooft calculation. The deviation of the 



'But the data for the upper states with the largest J are not reliable in some cases 
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trajectories from linearity may indicate a breaking of the 5*0(4) symmetry. Another possible 
explanation is that the effective value of N, in our formula appropriate for heavy mesons, 
is not iV — > oo (N = 100 used in Figs. 13 and 14), but much smaller. This may be due 
to coupling with break-up channels that effectively terminate the rotational and vibrational 
bands. The U(4) 3 0(4) gives the following dependence on N: 

M 2 (v) = Ml - 4(N + 1) A{v - v 2 /(N + 1)}, (16) 

with N ~ 20, we find the "vibrational" trajectories to be significantly bent, so as to actually 
describe the vibrational spectra of \l/ and T observed so far. More experimental work is 
needed to extend these studies of the excited vibrational states, to clarify the picture" (||, 
p. 910). 

We want to stress that few groups of authors simultaneously noted and discussed the 
nonlinearity of the experimental RT for J/\l/ and T families and its interplay with current 
theoretical quark models, namely Iachello ||, Semay |3(J, Stepanov 0, Dey and Sergeenko 

Em. 



In their series of papers Basdevant et al [381 employ semirelativistic hamiltonian and the 



HF method to describe spectra and RT of mesons and baryons. When authors have tried to 
fit simultaneously mesonic and baryonic spectra, they found that Regge slopes for mesons 
and baryons should be different. Authors describe experimental spectra of A, E, H and Q 
up to J = 19/2 with fair quality. 

Johnson et al [[39] used a semiclassical relativistic model for the orbital spectra of mesons, 



based on the assumption of a universal, flavor-independent linear confining interaction. They 
considered yrast lines for the families ir, 77, 7/, p, uj, K, K* , if, D, D* , D s , D*, and, in 
particular, noted that tt and K families do not lie on linear Chew-Frautschi plots. Johnson 
also found that the K* family and p family indeed have different slopes. 



Durand et al |4(J also obtained varying slopes, describing spin-averaged spectra of strange, 
charmed and bottom mesons in the Bethe-Salpeter approach. Their range of variation for the 
slopes is very similar to that of the ^.-expansion technique || , which is kind of a consistency 
check, because the two models started from different dynamical origins. 

Iachello jllj analyzed parity doubling phenomena in the meson and baryon sectors. He 



has found out that parity doubling (PD) definitely takes place at low momenta, maximally 
pronounced at J = 5/2, both experimentally and theoretically. Then, only at higher mo- 
menta linear RT appeared. "... One also observes that, at high angular momenta, parity 
doubling gradually disappears and linear Regge trajectories with slope identical to that of 
the meson spectra begin to appear..." (|41f, p. 38c, Figs. 10, 11). It is marvelous, that our 



model Jl0| practically reproduces independently Figs. 10, 11 from [41|, while employing 



NRQM with HF method. 

Sergeenko |J7|] derived a simple interpolating formula for the square of the quarkonium 
mass and an analytic expression for the RT in a whole region of both light and heavy quarko- 
nia on the basis of the consideration of two asymptotics for the QCD inspired interquark 
potential. He found effects of nonlinearity for the ^ and T families RT and also the effect 
of decreasing slope with the increasing quark mass. 

Very recently Filipponi et al || built a simple phenomenological model which attempted 
to describe analytically the quark-mass dependence in the mesonic RT for all flavors. These 
authors definitely got the gross feature of decreasing slope with increasing quark mass for 
RT in u, d, s, c, b - sector. The major drawback of the papers || is that the authors did 
not account for J-dependence of the slopes. In spite of this oversimplification, the trend of 



12 



their results roughly resembled that of our Table 3 for median values (ot')i over the whole J 
interval. 

Martin []IJ also considered parent RT for mesons and baryons in WKB-approach and has 
derived simple analytical formulas for flavor-dependent slopes, which are decreasing with 
increasing quark mass. Martin's formulas are very different from the Filipponi's results M, 
and he also noticed nonlinearity effects for RT in the resonance region. 

During the last few years group of T.Goldman obtained few new and very unusual results 
on RT f|3|]. Authors |4^] used an unquenched lattice potential and compute the spectrum 



of the bottomonium system. They demonstrate numerically that the effect of color screening 
is to produce termination of nonlinear hadronic parent and daughter's RT. These results are 
supported by recent lattice calculations, which observe string breaking [flBR . It would be 
interesting to study the effect of color screening on light quarkonia and baryon's RT. 

Burakovsky ||3|| constructed generalized string model for RT with the string tension 7 



depending on string point with broken reparametrizational symmetry. This model produces 
nonlinear Regge trajectories, which are in many cases analytical, such as square root, loga- 
rithmic and hyperbolic trajectories. 



L. D. Soloviev in Ref. [[mJ] suggested another string model with nonlinear RT admitting 
any growth rate for the trajectories. This power-law growth is limited by the exponent 3/2. 
Linear growth is typical for some special cases of quark interaction. The nonlinearity of the 



RT in the model J44j also results from some modification (dependence of the Lagrangian on 



scalar curvature) of the standard string action (|5|). In standard string models [23] the RT 
are nonlinear only for small J. 



There are very few papers considering RT for glueballs. Goldman et al [|!J used the 
"glueball dominance" picture of the mixing between qq mesons of different hidden flavors 
and established new glueball-meson mass relations, which predict the following glueball 
masses: M(0 ++ ) ~ 1.65 ± 0.05 GeV, M(l~-) ~ 3.2 ± 0.2 GeV, M(2~+) ~ 2.95 ± 0.15 GeV, 
M(3 ) ~ 2.8 ± 0.15 GeV. Their results are consistent with (quasi)linear RT for glueballs 
with a' G ~ 0.3 ±0.1 GeV" 2 . 



Yamada et al [pl| investigated RT of gluonia (2g) and hybrid mesons in the extended 
covariant oscillator quark model (COQM), and estimated glueball slope a' G ~ 0.26 GeV -2 . 
Authors predicted that the orbitally excited multiplets of gluonia with L = (1,2,3) lie in 
the mass region of (2.5 - 3.0 GeV, 3.2 - 3.6 GeV, 3.7 - 4.1 GeV). Their results are in fair 
agreement with lattice predictions. 

The relativistic quantum string model by L. D. Soloviev |TJ| describing reasonably well 
meson spectrum also predicts glueball RT with J PC = ++ — fo(1370) and fo (1500), which 
must lie on the following trajectories 



The results of our string and potential model fits and predictions for baryon and meson 
spectra and RT reveals distinctive feature — RT in many cases are nonlinear functions of J. 
This fundamental feature is in accord with analysis of pure experimental RT from PDG98 
(Table 5), and with predictions of dozens different quark models, reviewed in this paper. 
Regge trajectories for mesons and baryons are not straight and parallel lines in general 
in the current resonance region both experimentally and theoretically, but very often have 
appreciable curvature, which is flavor-dependent. 
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